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Abstract 

We show that a flowing dielectric medium with a linear response to an external electric field can be 
used to generate an analog geometry that has many of the formal properties of a Schwarzschild black 
hole for light rays, in spite of birefringence. The surface gravity of this analog black hole has a 
contribution that depends only on the dielectric properties of the fluid (in addition to the usual 
term dependent on the acceleration). This term may be give a hint to a new mechanism to increase 
the temperature of Hawking radiation. 

PACS numbers: 04.70.Bw, 04.70-s, 41.20.Jb 



I. INTRODUCTION 

In recent years, there has been a growing interest in 
models that mimic in the laboratory some features of 
gravitation 0|. The actual realization of these models 
relies on systems that are very different in nature: or- 
dinary nonviscous fluids, superfluids, flowing and non- 
flowing dielectrics, non-linear electromagnetism in vac- 
uum, and Bose-Einstein condensates (see for a com- 
plete list of references) . The basic feature shared by these 
systems is that the behaviour of the fluctuations around a 
background solution is governed by an "effective metric" . 
More precisely, the particles associated to the perturba- 
tions do not follow geodesies of the background spacetime 
but of a Lorentzian geometry described by the effective 
metric, which depends on the background solution. It is 
important to notice that only some kinematical aspects 
of general relativity can be imitated by this method 0, 
but not its dynamical features (see however 0>E|)- 

This analogy has permitted the simulation of several 
configurations of the gravitational field, such as worm- 
holes and closed space- like curves for photons 0,0. Par- 
ticular attention has been paid to analog black holes, be- 
cause these would emit Hawking radiation, as was shown 
first by Unruh in the case of dumb black holes 0- It is 
the prospect of observing this radiation (thus testing the 
hypothesis that the thermal emission is independent of 
the physics at arbitrarily short wavelengths 7]) that mo- 
tivates the quest for a realization of analog black holes in 
the laboratory. Let us emphasize that the actual obser- 
vation of the radiation is a difficult task from the point of 
view of the experiment, if only because of the extremely 
low temperatures involved. In the case of a quasi one- 
dimensional flow of a Bose-Einstein condensate for in- 
stance, the temperature of the radiation would be around 
70 nK, which is comparable but lower than the temper- 
ature needed form the condensate [a. 



In this paper we would like to explore the possibili- 
ties offered in the construction of analog black holes for 
photons by nonlinear flowing dielectrics p^ . Let us re- 
mark that, to our knowledge, all the articles that ana- 
lyze different aspects of black holes in dielectric fluids 
jl3l | published up to date are devoted to the case of a 
constant permittivity tensor. Here instead we would like 
to present a new static and spherically symmetric ana- 
log black hole, generated by a flowing isotropic dielectric 
that depends on an applied electric field. As we shall 
see, the radius of the horizon and the temperature of 
this black hole depend on three parameters (the zeroth 
order permittivity, the charge that generates the external 
field, and the linear susceptibility) instead of depending 
only on the zeroth order permittivity. Another feature 
of this black hole is that there is a new term in the sur- 
face gravity (and hence in the temperature of Hawking 
radiation), in addition to the usual term proportional to 
the acceleration of the fluid. This new term depends 
exclusively on the dielectric properties of the fluid, and 
it points to a new mechanism to get Hawking radiation 
with temperature higher than that reported up to date. 

We shall begin with the examination of the issue of 
photon propagation in a nonlinear medium in the next 
Section. We anticipate that photons experience birefrin- 
gence, i.e. there are two effective metrics associated to 
this medium. A photon "sees" one of these metrics or the 
other depending on its polarization. Although it may be 
argued that birefringence spoils the whole effective ge- 
ometry idea, we shall show in Section IIIII that photons 
see one and the same black hole, independently of their 
polarization. In Section Hvl we shall analyze in detail the 
case of a medium with a linear dependence in the exter- 
nal field. Some features of the motion of photons in this 
black hole are also exhibited in this section for the spe- 
cial case of a constant velocity flow, and it will be shown 
that there is a simple relation between the two effective 
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metrics. In Section [V] we calculate the effective surface 
gravity of these black holes for an arbitrary velocity pro- 
file, and display the appearance of the new term. We 
close in Section IvH with a discussion of the results. 



II. PROPAGATING MODES 

In this section we shall study the propagation of pho- 
tons with different polarizations in a nonlinear medium 
(see Ref.jl^ for details and notation). Let us define first 
the antisymmetric tensors F^i, and P^i, representing the 
electromagnetic field. They can be expressed in terms of 
the strengths {E, H) and the excitations {D, B) of the 
electric and magnetic fields as 

F^^y = v^^Ey - VyEf_, - rif_„^"'^VaBf3, 

where represents the 4-velocity of an arbitrary ob- 
server (which we will take later as comoving with the 
fluid). The Levi-Civita tensor introduced above is de- 
fined in such way that jy'^i^s _ Cartesian coordi- 
nates. Since the electric and magnetic fields arc spacelike 
vectors, the notation E°'Ea = -E^,H°'Ha = -H^ will 
be used. We will consider here media with properties 
determined only by the tensors and (i-e. me- 
dia with null magneto-electric tensor), which relate the 
electromagnetic excitations to the field strengths by the 
constitutive laws, 



D^^eJ{E,H)Ep, 



Br, 



nJ{E,H)Hp. (1) 



In order to get the effective metric, we shall use 
Hadamard's method The discontinuity of a func- 

tion J through a surface E will be represented by the 
symbol [J]s, and its definition is 

o — 

By taking the discontinuity of the field equations 
*F^"'.„ = and P^"'.y = 0, and assuming that 



and 



with fiQ = const., we get 



e(fc.e) - —{E.e){k.E) = 
E 



fj,o{k.h) — 



(2) 



(3) 



(4) 



(5) 



e(fc.w)e^ - ^E"er,ik.v)E^' + tj^""'^ k^Vah/s = (6) 
E 



Hoik.v)h'' - T^f'^'^^kyVc^ep = 



(7) 



where k^^ is the wave propagation vector, e' is the deriva- 
tive of e w.r.t. E, and 



[Ei_i,x] = kx, 



[Hf_,^x] = kx 



Note in particular that Eqn. Ql shows that the vectors k^ 
and e'^ are not always orthogonal, as would be the case 
if e' was zero. 

Substituting Eqn.Q in ©, we get 

Z^'^g^ ^ (8) 

where the matrix Z is given by 

= [k^ + {k.vf{fioe-l)]v''^ 

-^^^.(k.v)^E^'E'^ + {v.k){v''k'^ + fc^w'^) 
E 

- [e^io{k.v) + v>'v^ ~ k^k^. (9) 

Non-trivial solutions of Eqn.|(Sl) can be found only for 
cases in which det | Z'^^ | = ( this condition is a gener- 
alization of the well-known Fresnel equation pH'). 

Eqn.|(Sl) can be solved by expanding e^, as a linear com- 
bination of the four linearly independent vectors Vi,, Ei,, 
ki, and r/apfj.iyV^E^k'^ (the particular case in which the 
vectors v^, E^, and k^, are coplanar will be examined be- 
low). That is. 



Zy = aEy + Pr]aXf,yv"E^k'' + -fk^ + Sv^. 



(10) 



Notice that taking the discontinuity of E^.^ we can show 
that (e.v) = 0. This restriction imposes a relation be- 
tween the coefficients of Ean. lfTT 



S ~ —"f{k.v) 

With the expression given in Ean. llO|l . Eqn. reads 

1 



a[k^ -{l-fio (e Ey){k.vf]-j 



ti^ik.vY- e'E'^kr, 



= 



aS^fc^ + 7(1 - HQe){k.vf + 6{k.v) = 
a{k.v)E^'k^ + -i{k.v)k^ + 5[k^ + fxoe (k.v)^] = 
(3[k^-{l-fioe){k.vf] =0 

The solution of this system results in the following dis- 
persion relations: 

kl = ik.v)^^-^lo{eEy] + ^e' E''E'^ko,kp,{ll) 
kl = [1 - fioe{E)]{k.v)\ (12) 
They correspond to the propagation modes 

e; - p- {^ioe{k.vyE, + E''kr,[k,-{k.v)v,]} 



et = P+ VcXt^uv'^E^k^', 



(14) 



where and p+ are arbitrary constants. The labels 
and "— " refer to the ordinary and extraordinary rays, re- 
spectively. Eqns. (|ll|l and p2|) govern the propagation of 
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photons in the medium characterized by ^ = /iq =const., 
and e = e{E). They can be rewritten as g'^k^^k^ — 0, 
where we have defined the effective geometries 

.9^ = V^"" - [1 - Aio (e E)'] v^v'^ - ^ e' E^E'^ilb) 
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(+) 



(16) 



The metric given by Ean. (|15|l was derived in [T^ l. while 
the second metric very much resembles the metric de- 
rived by Gordon The difference is that in the case 
under consideration, e is a function of the modulus of 
the external electric field, while Gordon worked with a 
constant permeability. 

Let us discuss now a particular instance in which the 
vectors used as a basis in Eqn. (fTn|) are not linearly inde- 
pendent. If we assume that 



(17) 



then, by Ean. (|10|l . the vectors e^, k^, and are copla- 
nar. In this case, the basis chosen in Ean. HlU|) is not 
appropriate. Notice however that if we assume that 
is a combination of vectors that are perpendicular to fc^, 
so that (e.fc) = 0, Eqn.lO impHes that (E.e) = 0. The 
converse is also true: if {E.e) = 0, then, from Ean. (|17|l . 
(fc.e) — 0. For this particular case, in which is per- 
pendicular to u'^, fc'^ (and consequently to -E^), Eqns.© 
and lO imply that 

[fc2 + (fc.z;)2(^oe-l)] e^ = 

We see then that in the case in which E'^ — ak^ + bv^, 
Fresnel's equation determines that the polarization of the 
photons is perpendicular to the direction of propagation 
and to the velocity of the fluid. Moreover, the motion 
of these photons is governed by the metric g'i^ . For in- 
stance, if the electric field, the velocity of the fluid, and 
the direction of propagation are all radial, then the polar- 
ization is in the plane perpendicular to the propagation, 
and the two polarization modes feel the same geometry. 



III. THE ANALOG BLACK HOLE 

We shall show in this section that the system described 
by the effective metrics given by Eans. H15|l - H16|) can be 
used to produce an analog black hole. It will be conve- 
nient to rewrite at this point the inverse of the effective 
metric given by Ean. ljlSIl using a different notation: 



■ - ' 1 + ^ 
where we have defined the quantities 
1 



lt.lu, (18) 



/ 



E, 

e ' - E- 



Note that e — e{E). We have introduced here the velocity 
of light c, which was set to 1 before. 



Taking a Minkowskian background in spherical coor- 
dinates, and 

= (vo,fi,0,0), E^ = (£;o,Ei,0,0), (19) 

we get for the effective metric described by Ean. (|18|l . 



(-) ^"^1 n ^^ , ; ; 



(20) 
(21) 
(22) 



and (?22 "* ^i^'i 533 "* Minkowski spacetime. The vec- 

tors and satisfy the constraints 



2 2 2 

t;„ - = C , 

vqIo - vih ^ 0. 



(23) 
(24) 
(25) 



This system of equations can be solved in terms of vi, 
and the result is 



^1; 



7)2 

;2 _ ^1 
" C2 ' 



/2 - 

n — 



c2 + 



(26) 



(27) 



Now we can rewrite the metric in terms of /3 = vi/c 
[23 |. The explicit expression for the metric coefficients is: 



(_) _ l-/32(c2^oe-l) 
c^M^ + e'E) ' 



5oo 



9o 



1 - c^^pe 
c2^o(e + e'E)' 



(28) 
(29) 
(30) 



From Ean. (|28fl it is easily seen that, depending on the 
function e(E), this metric has a horizon at r — r^^ given 
by the condition goo{rfi) = or, equivalently. 



(31) 



The metric given above resembles the form of 
Schwarzschild's solution in Painleve-GuUstrand coordi- 
nates [mil: 



2GM\ 



dt^ ± 2 



2GM 



dr dt- dr^ -r^dn\ 
(32) 
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With the coordinate transformation 



j2GM/r , 
dtp — ats T TT7TTT- dr, 



1 



2GM 



(33) 



the hne element given in Ean. H32|l can be written in 
Schwarzschild's coordinates. The "+" sign covers the 
future horizon and the black hole singularity. 

The effective metric given by Eans. (|28f) - (|3UII looks like 
the metric in Ean. l|32() 25] . In fact, it can be written in 
Schwarzschild's coordinates, with the coordinate change 



dt 



PG 



dts 



.901 W 
goo{r)' 



dr. 



(34) 



Using this transformation with the metric coefficients 
given in Eans. H28|l and H29|) . we get the expression of 
g[i^ in Schwarzschild coordinates: 



fii 



eiE) 



il-p^c^fioeiE)-l]ME) + eiEyE)- 



(35) 



Note that ^ is zero in the new coordinate system, while 

(7oo '' is still given by Eqn. H28|l . Consequently, the posi- 
tion of the horizon docs not change, and is still given by 
Eqn.(ETJ- 

Working in Painlcvc-GuUstrand coordinates, we have 
shown that the metric for the "— " polarization describes 
a Schwarzschild black hole if Eqn.(j2U has a solution. Af- 
terwards we have rewritten the "— " metric in the more 
familiar Schwarzschild coordinates. Let us consider now 
photons with the other polarization. They "see" the met- 
ric given by Ean. (|16|l . whose inverse is given by: 



5^1/ -Vt^i' c c ' 



1 



(36) 



Using this equation and Eans. (|26|l and H27() it is straight- 
forward to show that 



yoo 



V c Moe(£^)y 



1 - 



c^f^oe{E) 



(37) 



(38) 



1 



l + /32(l-c2/ioe(£^))' 



(41) 



It is important to stress then that the horizon is located 
at rh given by Eqn. L'U\I for photons with any polarization. 
Moreover, the motion of the photons in both geometries 
will be qualitatively the same, as we shall show below. 



IV. AN EXAMPLE 

We have not specified up to now the functions e{E) and 
E{r) that determine the dependence of the coefficients of 
the effective metrics with the coordinate r. From now on 
we assume a linear e{E) j26|., 

e(E)^eo(x + X^^^Eir)), (42) 

X^^^ ■ The nontrivial Maxwell's equation 



with X = 1 
then reads 



(V^e(r)i^°i),^=0, 



(43) 



where 7 is the determinant of the flat background metric. 
Taking into account that 



p2 



(44) 



we get as a solution of Ean. (|43|l for a point source in a 
flat background in spherical coordinates 



F 



01 



-X± Vx" + 4x(2)g/£or2 
2cx(2) 



(45) 



Let us consider a particular combination of parameters: 
X*-^^ > 0, Q > and the sign in front of the square 
root in F'^^ , in such a way that i? > for all r. To get 
more manageable expressions for the metric, it is conve- 
nient to define the function a{r): 



E{r) 



X 



2x(2 



cr(r) 



(46) 



where 



9ii 



1 



(39) 



This metric also corresponds to a Schwarzschild black 
hole, for some e{E) and /3. Comparing Eans. (|28|l and 
(|37|l we see that the horizon of both analog black holes 
is located at rh, given by Eqn.|j3lJ- 

By means of the coordinate change defined by 
Ean. (|34|l . we can write this metric in Schwarzschild's co- 
ordinates. The relevant coefficients are given by 



(+) _ i+^^O-cWGB)) 



(40) 



and 



cr(r) 



-1 + -y/r^ +q 
r 



4x( 



2)/ 



In terms of cr, the metrics take the form 



dr' 



2-/32 [x(a(r)+2)-2] 
2x (l + fT(r)) 

2 + a(r) 

[2 - /32 (x ia{r) + 2) - 2)] (1 + a{r)) 



(47) 



(48) 



(49) 



r^dn^ 



5 



2-/3" [xHr)+2)-2] 
x(2 + <T(r)) 
2 

2 + /32 [2-x(a(r)+2)] 



dr' - (50) 



Notice that the {t, r) sectors of these metrics are related 
by the following expression: 



where the conformal factor $ is given by: 

1 + cr(r) 



(51) 



$ = 2 



2 + a{r) 



We shall study next some features of the effective black 
hole metrics. We would like to remark that up to this 
point, the velocity of the fluid vi is completely arbitrary; 
it can even be a function of the coordinate r. We shall 
assume in the following that vi is a constant. This as- 
sumption, which will be lifted in Sect0 may seem rather 
restrictive but it helps to display the main features of the 
effective metrics in an easy way. 

To study the motion of the photons in these geome- 
tries, we can use the technique of the effective potential. 
Standard manipulations (see for instance [2^) show that 
in the case of a static and spherically symmetric metric, 
the effective potential is given by 



V{r) 



1 



1 



gm{r)gii{r)J r^gii{r) 



(52) 



where e is the energy and L the angular momentum of 
the photon. 

In terms of cr(r), and of the impact parameter = 
L^/e^, the "small" effective potential v{r) = V{r)/e-^ for 
the metric Eqn.jSDI) in Schwarzschild coordinates can be 
written as follows: 



(r) = l 



2(l+a(r))2 6^ (2-/3V(r))(l + a(r)) 
2 -I- cr(r) r2 



r2 2 + (T(r) 

(53) 

A short calculation shows that v'^ Ms a monotonically 
decreasing function of p. Consequently, we shall choose 
a convenient value of it, for the sake of illustrating the 
features of the effective potential. Figures (Q) and 
show the plots of the potential for the (— ) metric for 
several values of the relevant parameters. 

The effective potential for the Gordon-like metric can 
be obtained in the same way. From Eqns.jS^J and 
we get 



/+\r) = 1 - 



cr(r) 



— [2-/3MO] (54) 



The plots in Figures Q and show the dependence of 
v(+)(^r) on the different parameters. 

We see from these plots that, in the case of a constant 
flux velocity, the shape of the effective potential for both 
metrics qualitatively agrees with that for photons moving 
on the geometry of a Schwarzschild black hole (see for 
instance Ref.|2^, pag. 143). 
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FIG. 1: Plot of the effective potential \r) for g = 1, 6 : 
1,3,5 (starting from the lowest curve), and /3 — 0.5. 
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FIG. 2: Plot of the effective potential \r) for 6 = 3, and 
q — 1, 3, 5 (starting from the lowest curve), and /3 — 0.5. 



V. SURFACE GRAVITY AND TEMPERATURE 

Let us now go back to the more general case oi /3 — 
f3{r), and calculate the "surface gravity" of our analog 
black hole. We present first the results for the constant 
permittivity case. By setting e'{E) = in the metrics 
Eans. H15|l and (|16|l . we regain the example of constant 
index of refraction studied for instance in 13] . It is easy 
to show from Ean. H31|) that the horizon of the black hole 
in this case is given by 



1 



X-1 



(55) 
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FIG. 3: Plot of the effective potential for the Gordon-like 
metric, for q = 1, b = 1,3, 5 (starting from the lowest curve), 
and p = 0.5. 




FIG. 4: Plot of the effective potential for the Gordon-like 
metric, for 6 = 3, g = 1, 3, 5. 



for 0^{rh) ^ 1. This acceleration is a quantity that must 
be determined solving the equations of motion of the fluid 

[23. 

Going back the the more general case of a linear per- 
mittivity, described by the metrics given by Eans. H50() 
and (|5T|l . and considering that (3{rh) ^ 1, the radius of 
the horizon is 1281: 



qx 



(58) 



Using the expressions given above, the result for the sur- 
face gravity of the "— " black hole for P{rh) ^ 1 is 



1 



(59) 



This equation differs from the surface gravity of the case 
of constant permittivity (Ean. H57|) ) by the presence of a 
new term that does not depend on the acceleration of the 
fluid. To see where this new term comes from, we can 
take the derivative of the metric coefficient 1700 ) given 
by Ean. (|28|l . with respect to the coordinate r (see the 
definition of k, Ean. (|56|l . Taking into account Ean. (|42|l . 
the resuh is 



500,r 



1313' (ejE) + ejEYE) (1 - c^fipejE)) 
cVo ieiE)+eiEyEr 



dt{E) 



In this expression, the first term is the generalization of 
the term corresponding to the case e = const, (compare 
with Ean.H57|l). which mixes the acceleration of the fluid 
with its dielectric properties. On the other hand, the 
second term, which is the new term displayed in Ean. (|59|l . 
is related only to the dielectric properties of the fluid. 

It is easy to show that the these results also apply 
to the black hole described by the Gordon-like metric. 
This is not surprising though, because of the conformal 
relation between the two metrics, given by Ean. (|51|l (22l |. 



The "surface gravity" of a spherically symmetric analog 
black hole in Schwarszchild coordinates is given by [l7| 



K = — lim 



.900,: 



2 ''-^"^ ^/\gll\ goo 



(56) 



For the metrics Eans. l|15|l and (|16|l with e = eox and 
given by Ean. (|55|l . the analog surface gravity is 



2 Vx ^' ' 



(57) 



In this expression we can see the influence of the dielectric 
properties of the fluid (through the constant x) and also 
of its dynamics through the physical acceleration in the 
radial direction, given by 



2 ^ ' ii-h 



VI. DISCUSSION 

We have shown that a flowing inhomogeneous dielec- 
tric that depends linearly on an external electric fleld 
generates one effective metric for each polarization state 
of photons. A particular conflguration of the fluid ( i.e. 
pure radial flow) plus an external electric fleld was shown 
to be an analog black hole, with a radius that depends 
on the function (-(E). The features of the black hole de- 
pend on the charge that generates the electric fleld, on 
the properties of the dielectric, and on the velocity pro- 
flle. Note in particular that in the absence of flux (/3 = 0) 
the metrics do not describe a black hole. 

The existence of two metrics reflect the birefringent 
properties of the medium. Although some claims have 
been made that birefringent materials spoil from the be- 
ginning the idea of an effective metric, we have shown 
here that for a special conflguration, even if two metrics 
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are present, photons with different polarizations experi- 
ence the same horizon. Moreover, as seen from the plots 
of the effective potential, the motion of these photons 
in the medium will depend on their polarization, but is 
qualitatively the same for both types of photons. 

At first sight it may seem that by choosing an appro- 
priate material and a convenient value of the charge we 
could obtain a high value of the temperature of the radi- 
ation, given by |29j, 

T = — ^ K « 4 X 10"^^ K KsVm, (61) 
2'KkBC 

with K given in Ean. (|59|l . However, the equation for the 
surface gravity can be rewritten as 

»-'(^^')L■ 

We see then that, because (3{rh) ^ 1, the new term 



appearing in k is bound to be very small, but of the 
same order of magnitude of the acceleration term. To 
summarize, we have shown that a new term emerges in 
the surface gravity which is produced by the properties 
of the material and is comparable in magnitude to the 
acceleration term. This result indicates that it may be 
worth to study if some media with nonlinear dependence 
on an external electric field are appropriate to generate 
analog black holes with Hawking radiation with a higher 
temperature than that produced up to date with other 
systems. 
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